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ABSTRACT 



The stellar wind around compact object in luminous wind-accreting high-mass X-ray binaries is expected to be strongly ionized with the 
X-rays coming from the compact object. The stellar wind of hot stars is mostly driven by light absorption in lines of heavier elements, 
and X-ray photo-ionization significantly reduces the radiative force within the so-called Stromgren region, leading to wind stagnation 
around the compact object. In close binaries like Vela X- 1 this effect might alter the wind structure throughout the system. Using the 
spectral data from the Monitor of All-sky X-ray Image (MAXI), we studied the observed dependence of the photoelectric absorption as 
a function of the orbital phase in Vela X- 1 and found that it is inconsistent with expectations for a spherically symmetric smooth wind. 
Taking into account previous investigations, we developed a simple model for wind structure with a stream-like photoionization wake 
region of slower and denser wind trailing the neutron star, which is responsible for the observed absorption curve. 
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Vela X-1 is a persistently active high-mass X-ray binary system 
consisting of a massive neutron star (1.88M Q , Quaintrell et al. 
|2003[ ) and the B0.5Ib-type supergi ant HP 77581 with a mass 
of ~ 23M G and radius of ~ 30R o ( |van Kerkwijk et al |[l995l 



The hot primary loses mass at a rate of ~ 10 °M Q yr 1 (|Nagase 



|et al.|1986| > through a powerful stellar wind, which is accreted 
onto the neutron star located ~ 53R Q away, and powers the X- 
ray pulsar. The pulsar is eclipsed by the primary during every 
orbital cycle of ~ 8.964 d ( |van Kerkwijk et al.|1995). T he spin 
period of the neutron star, P s ~ 283.5 s (Rappapo rt|1975| l, has re- 
mained almost constant since this discovery. The X-ray spectrum 
of Vela X-1 is typical for accreting pulsars, and is well described 
by a cutoff power law (Nagase et al. 1986) with cyclotron reso- 
nance scattering feat ures at ~ 25 ([ Makishima & Mihara 1992) 
and ~ 50 - 55 keV ( |Kendziorra et al.||1992[|Maldshima et al.| 
1999). Similar to other wind-accreting systems, the X-ray flux in 
Vela X-1 is strongly variable and may change by orders of mag- 
nitude within minutes ( Kreykenbohm et al. 2008 1 Doroshenko 
et al. 201 1 ). The variability is usually associated with the clumpy 



nature of the captured radiatively driven wind (Oskinova et al. 
|2007|[Furst et al.|2010| >. The observed average X-ray flux implies 
a luminosity of ~ 4 X 10 36 erg s _1 (assuming a distance of ~ 2kpc, 



Nagase 1989), which agrees with the estimate for Bondy-Hoyle 
accretion from smooth spherically symm etric wind with an ob - 
served terminal velocity of ~ 1 100 km s ( jWatanabe et al. 2006), 
and a mass-loss rate of ~ 1O~ 6 M yr _1 . The wind in Vela X-1, 
however, is not expected to be symmetric, and the X-ray emission 
from the pulsar is supposed to ionize the wind almost down to the 
surface of the primary, leading to wind stagnation and formation 
of the photoionization wake trailing the neutron star ( |Blondin| 



et al. 1 1 990 Feldmeier et al. 1996). Similar studies were also per- 



formed for other sources (Manous akis et al.|20 12). More recently, 
detailed calculations of wind acceleration in the presence of an X- 
ray source have been carried out by |Krticka et al. ( 2012| i with ap- 
plication to Vela X-1 for the stationary case. The authors found 
that the wind velocity field is significantly affected by the feed- 
back from the accreting neutron star and anticipated an almost 
complete stagnation of the wind within the "photoionization- 
bubble" around the neutron star. 

Several studies point to the existence of such a feature in 
Vela X- 1 . Nagase et al. ( 1986 1 observed strong asymmetry be- 
tween the ingress and egress of the eclipse, accompained by in- 
crease in energy of Fe-K line and the absorption edge, which 



Feldmeier et al. ( 1996 1 interpreted as enhanced scattering of X- 



rays in the denser region, trailing the neutron star. Based on the 
large observed size of the scattering region, these authors ruled 
out the accretion wake in favor of the photoionization wake as 
the scattering site. 

Based on the orbital modulation of the UV- resonance lines 
in several systems, including Vela X-1, van Loon et al. ( 2001| l 
concluded that a complex wind structure with an ionized region 
around, and a photoionization wake behind the neutron star is 
required to explain the high-resolution UV spectra of the source. 

Goldstein et"aL] ( |2004| i and |Watanabe et al.| ( |2006) used high- 



resolution Chandra X-ray spectral data to constrain the ioniza- 
tion structure and the geometrical distribution of material in the 
system. Intensities of the X-ray emission lines of highly ionized 
H- and He-like ions were found to be compatible with the smooth 
wind model ( |Castor et al.|1975) . However, most of the X-ray re- 
combination lines were found to come from between the neutron 
star and the primary, suggesting an hanced density in this region, 
most likely due to photoionization-driven reduction in radiative 



V. Doroshenko et al.: Footprints in the wind of Vela X-l traced with MAXI. 



force. In fact, the velocity shifts of the highly ionized recombi- 
nation emission lines were found to be smaller in magnitude by 
several hundred kilometers per second than the predictions for 
the smooth wind,which directly supports the presence of the pho- 
toionization wake. 

Performing a detailed study of the global wind structure is 
complicated, however, by the dramatic variability of both the 
wind, and the X-ray feedback from the pulsar. Studying this vari- 
ability provides insight on details of wind acceleration and interac- 
tion with the neutron star and is best addressed with high-quality 
dedicated observations such as those presented for Vela X-l by 
Watanabe et al. ( 2006) or[Nagase et al. ( 1986 1, who performed 



detailed time-resolved spectral analysis. On the other hand, to to 
trace the persistent wind structure, not only long observations at 
different orbital phases, but also observations at multiple orbital 
cycles are important, which might not be feasible for dedicated 
observations. On the other hand, because it is one of the brightest 
persistent X-ray sources in the sky, Vela X-l is an ideal target for 
monitor-type X-ray instruments such as MAXI (Monitor of All- 
sky X-ray Image), which can complement existing high-quality 
dedicated observations. In this work we used MAXI data to study 
the average orbital phase dependence of the X-ray spectrum in 2- 
30 keV in Vela X- 1 with the aim of constraining the changes in 
photoelectric absorption induced by the photo-ionization-driven 
inhomogeneties in wind structure. We found that the observed 
changes in the absorption column are incompatible with the spher- 
ically symmetric smooth wind, and argue that the photoioniza- 
tion wake trailing the neutron star may indeed be responsible for 
the enhanced absorption. 

2. Observations and data analysis 

MAXI is the first astronomical mission operated on the 
International Space Station (ISS) (Matsuoka et al.|2009| . It was 
installed on the Japanese Experiment Module — Exposed Facility 
on July 16, 2009, and has been designed to surpass all operating 
X-ray monitors both in sensitivity and energy resolution. Two slit- 
slat X-ray cameras scan the entire sky every 92 min as the ISS 
follows its orbit. They operate in 0.5-12 keV (the Solid-state Slit 
Camera (SSC) |Tomida et al.|201 T) , and 2-30 keV (the Gas Slit 
Camera (GSC), Mih ara et al. 201 1 ), which we used here. The cam- 
era consists of twelve Xe-gas proportional counters with slit-and- 
slat collimators limiting the field of view (FOV) to 1 .5° x 160°, 
which scans the sky as the ISS orbits. The counters employ re- 
sistive carbon-wire anodes to acquire one-dimensional position 
sensitivity, and enable position reconstruction for every source 
in the FOV to about 30'. The large detector area of 5350cm 2 al- 
lows many brighter sources to be detected within one ~ lm scan. 
Typical daily exposures reach 3000 cm 2 s, which enables spectral 
studies when multiple scans are combined. The relatively low un- 
certainty in the energy calibration of about 3%, the fair energy 
resolution ~ 18 % (at 5.9 keV), and the robust background rejec- 
tion via anti-coincidence with outer veto cells additionally aid 
spectral capabilities. The camera performance is summarized in 
Sugizaki et~alJ ( |2T)TT] ). 

The inspection of public Vela X- 1 light curves provided by 
the MAXI teanrlfor every source detected by the GSC in several 
energy ranges suggests a significant dependence of the source 
spectrum on orbital phase. The comparison of the mission-long 
soft and hard light-curves folded with orbital period presented in 
Fig.[T]reveals that most changes occur in the soft 2-4 keV band, 
while other bands are less affected, which could be naturally 
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Fig. 1. Vela X-l orbital light curves in 2-4 (blue circles), 4-10 
(green triangles), and 10-20keV (red squares) energy ranges as- 
suming orbital parameters reported by Kreykenbohm et al. (2008 ), 
and folded from public MAXI light curves relative to the perias- 
tron. The arrow indicates the mid-eclipse phase. 



explained with the strong changes in absorption column along the 
orbit. It is interesting also to note the spike in hard X-ray flux at 
phase ~ 0.4, which could potentialy provide additional constrains 
on the wind structure. 

To go beyond studying the hardness ratios and to verify if 
their change is related to the photoelectric absorption, we per- 
formed an orbital -phase-resolved spectral analysis of MAXI GSC 
data obtained during period MJD 55141-56200 using the standard 
MAXI pipeline (Sugizaki et a l.|2011||Nakahira et al.|2012| > modi- 
fied to allow arbitrary grouping of the individual scans. The indi- 
vidual scan duration is just about one minute, i.e., much shorter 
than the phase bin of ~ 1 day we are interested in, so grouping 
of the scans is a robust way to obtain the orbital-phase-resolved 
spectra. Using the ephemeris provided by Kreykenbohm et al. 
(2008 ), we calculated a list of start and end times for each phase 
bin, and extracted the spectra by grouping the individual scans ac- 
cordingly. The effective exposure is uniformly distributed along 
the orbit and varies between 15-20 ks per phase bin. We rebined 
all extracted spectra to contain at least 200 photons per energy 
bin, and fit them in the 2-20 keV energy range using the comp- 
tonization model by Sunyaev & Titarchuk ( 1980 1 (CompST in 
XSPEC). We kept the temperature and optical depth the same 
for all spectra, and only allowed the flux and the absorption col- 
umn to vary independently. This resulted in a statistically accept- 
able combined fit with;^ 2 ^ = 1.13 (x - 941.2 for 831 degrees 
of freedom) with the best-fit values for temperature and optical 
depth kT - 6.4(5) keV, r = 13.3(7), and the absorbtion column in 
3.3(5) - 25(4) x 10 22 atoms cm 2 range. Estimating the continuum 
parameters independently does not improve the fit significantly. 

Keeping the continuum parameters common for all phase 
bins is also justified because the observed shape of the spectrum 
does not change above lOkeV (within available statistics). We 
also confirmed that the data point distribution on a color-color 
plot is consistent with them being one-dimensional, controlled 
only by the absorption. The representative spectra and the orbital 
dependence of free parameters are presented in Fig. [2] The flux 
plotted in the figure is unabsorbed source flux in the 2-20 keV 
range derived from the spectral fit, and uncertainties are given at 
lcr confidence level. 
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Fig. 2. Representative unfolded spectra for the most (red) and 
least (black) absorbed phase bins (top panel), and best-fit results 
for the absorption column and flux for all spectra (bottom panel). 
Absorption column and flux estimates for the smooth wind (black) 
and the orbital light curve (gray) are plotted in the bottom panel. 
Model flux mostly changes because of change in the relative 
velocity of the wind and the neutron star. 



3. Discussion 

First of all we would like to emphasize that the observed or- 
bital dependence of the absorption column is asymmetric with 
respect to the conjunction, and thus cannot originate in a spheri- 
cally symmetric wind of the primary. The asymmetry can only 
be caused by interaction of wind with the neutron star. The ob- 
served dependence closely resembles the results of hydrodynam- 
ical simulations carried out by |Blondin et al.| ( |1990| l that ac- 
count for a photo-ionization wake (see for instance Figs. 8 and 
1 1 in their work), although the observed asymmetry is somewhat 
stronger in Vela X- 1 . A similar orbital dependence of the ab- 
sorption column was observed in IGR J17252-3616 with XMM- 
Newton (Manousakis & Walter 201 1 ), where the authors also at- 
tributed the enhanced absorption at later orbital phases to a photo- 



ionization wake. For the same source, Manousakis e t al.| ( [20T2] l 
performed detailed hydrodynamical simulations that confirmed 
their initial assumption that the peculiarly slow wind and the 
photo-ionization wake trailing the neutron star are responsible 
for the observed orbital dependence of the absorption column. 

Indeed, comparing the observations with the results of de- 
tailed hydrodynamical simulations is essential for understanding 
the effects the accreting neutron star imposes on the wind struc- 
ture. However, one needs to bear in mind the extreme complex- 
ity of these simulations and associated limitations. For instance, 
Manousakis et al. (2012) estimated the mass of the neutron star 
based on its impact on the simulated wind structure, which may 
be affected by the assumed simplified prescription for the radia- 
tive force that drives the wind (which is assumed to switch off 



within the Stromgren sphere). On the other hand, Krticka et al. 
(2012) showed that the radiative force field is much more compli- 
cated than that and has to be calculated self-consistently with the 
wind structure. Performing such calculations is beyond the scope 
of the present work, and taking into account that the results of a 
simplified approach presented in Manousakis & Walter (2011 1 
and the more detailed simulations presented in Manousakis et a . 
(2012) agree reasonably well, we stick to the simpler prescription 
just to illustrate some important points suggested directly by the 
observations. 

As a first step, it is interesting to compare the observed ab- 
sorption column and flux with estimates from a standard smooth 
radiatively driven wind model I Castor et al.|1975| > with the radial 
wind velocity profile, 



u(r) = K + (v a 



«>(.-* 



where vo ~ 10 km s is the sound velocity at the surface of the 
optical star with radius R s — 30R Q , v^ — 1100 km s ~ and/? ~ 1 
(Watanabe et al. (2006) and references therein). For the geom- 
etry of the system we assumed the orbital parameters reported 
by Kreykenbohm et al. ( 2008) , that the system is observed edge 
on ( |van Kerkwijket al. ( 1995 i esti mate ; > 75°), and a dopted a 
mass-loss rate of 2 x 10 M o yr _1 ( [Watanabe et al.|2 006). The 
density at any distance may then be calculated from the continuity 
equation. This allows one also to calculate the X-ray flux assum- 
ing the Bondi-Hoyle spherical accretion and the distance to the 
source of 2kpc (Nagase et al.|1986) . The equivalent absorption 
column was calculated by integrating the density along the line of 



sight and assuming neutral absorber cross-sections (Morrison & 
McCamm on|1983 1 with solar abudances. We neglected interstel- 
lar absorption in the direction of the source, which is e stimated to 
be < 5 x 1 21 atoms cirT 2 (Dickey & Lockman 1990), i.e., about 
an order of magnitude lower than the lowest observed value. The 
results are presented in Fig. [2] 

As previously mentioned, this simple model cannot reproduce 
the observed asymmetric shape of the absorption curve along 
the entire orbit. On the other hand, if the observed asymmetry is 
caused by some structure trailing the neutron star, the agreement 
between the observed absorption column and the model for the 
first three points after the eclipse when this structure is behind 
the neutron star is remarkable. The observed average flux agrees 
well with our simple estimate, although its orbital dependence 
is not reproduced as well. It is important to emphasize that no 
parameters were adjusted to obtain the curves in Fig [2] which 
were instead calculated using existing independent estimates for 
the parameters of the wind and the binary system. 

Another important point is that both the absorption column 
and the flux depart from the model predictions, and that the 
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observations start to depart from the model before the conjunction 
(phase 0.5). This means that the neutron star is not simply trailed 
by the region of slower, denser wind, but is submerged within it or 
at least resides at its boundary. This suggests a possible additional 
gradient in wind velocity in the vicinity of the neutron star, which 
might have profound effects on the amount of angular momentum 
captured by the neutron star, and therefore on its spin evolution. 
Detailed simulations with proper treatment of the radiative force 
and accurate wind velocity calculations are required to estimate 
their magnitude, however. 

To better illustrate these points, and to qualitatively explain 
the observed changes in the absorption column along the orbit, we 
extended a simple smooth-wind model by adding a bent stream 
of slower, denser wind trailing the neutron star similar to the one 



emerging from simulations by Blondin et al. ( 1990), Feldmeier 
et al.| ( 1996 1, and Manousakis et al. (2012|). Even these quite 



detailed simulations probably do not provide a realistic picture, 
since the accelerating radiative force is treated very crudely. For 
our "toy" model we just assumed that the expected stream-like 
structure is described as a wind sector with a certain opening an- 
gle (~ 40°), the same mass-loss rate and velocity profile as the 
ambient wind. The wind velocity within the stream was, how- 
ever, scaled down by a constant factor, thus increasing the density 
in the stream with respect to the ambient wind by about a fac- 
tor of ten, which is calculated from the continuity equation. The 
stream then bends away from the neutron star and phase lag lin- 
early increases with radius (oc O.lr), which corresponds to rigid 
rotation (Feldmeier et al.|1996 1. In reality, the wind velocity is 
certainly more complex, particularly between the companions 
(which implies, for instance, that the density wake will be offset 
from the primary). However, the line of sight does not cross this 
region, and most of the absorbing material resides in the outer 
wind, which has a simpler structure. The absorption curve calcu- 
lated using our "toy" model is qualitatively consistent with obser- 
vations as presented in Fig. [3] and the corresponding animation^] 
The result is also quite similar to the one reported by Manousakis 
& Walter| pOTT) for IGR J17252-3616, which is not surpris- 



ing because they used similar assumptions. We emphasize that 
in our case the model has mostly illustrative purpose, and only 



in conjunction with more detailed simulations (i.e., Manousakis 
et al.|2012 1 it is possible to conclude that the absrobtion in the 



dense stagnated wind region in the photoionization wake is most 
likely responsible for the observed absorption orbital profile in 
Vela X-l as well. We did not attempt to model the orbital flux 
profile here, since the accretion rate mostly depends on the rela- 
tive velocity of the wind and the neutron star, which cannot be 
estimated without detailed calculations of the wind-driving force. 
Taking into account the reasonable agreement of the observed 
flux and estimated luminosity for unperturbed wind (Fig. El, it is 
clear, however, that the wind velocity in the vicinity of the neu- 
tron star must not drascically decrease, i.e., it most likely resides 
at the edge of the photoionization wake. 

4. Conclusions 

We have investigated the global stellar wind structure in Vela X-l 
exploiting the continuous monitoring of the source with MAXI 
and its admirable spectral capabilities. We measured the average 
orbital dependence of the absorption column in the system and 
of the intrinsic source X-ray flux in Vela X- 1 , and found that 
both depart significantly from the expectations for a smooth-wind 
model. We attributed these discrepancies to a photo-ionization 
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Fig. 3. Top panel: Sketch of the system with the stream trailing 
the neutron star as it orbits HD 77581 (top panel); the solid line 
indicates the wind density distribution along the line of sight 
(indicated with an arrow). Bottom panel: orbital absorption curve 
as observed with MAXI (crosses, uncertainties at lcr confidence 
level), and model curves calculated for the smooth-wind model 
(dashed line, see also Fig.[2ji, and for a model with stream trailing 
the neutron star (solid line). We arbitrary shifted the curve for the 
smooth model to enhance the plot clarity. 



wake, which is expected to manifest as a denser stream-like re- 
gion trailing the neutron star. We constructed a simple model to 
qualitatively describe the impact of this stream on the orbital de- 
pendence of the absorption column, and found that the neutron 
star must reside either within it or at the leading edge. We found 
that this simple mode agrees quite well with observations, which 
provide important constraints for future more detailed wind simu- 
lations. From an observational point of view, our result also poses 
a direct science case for and highlights the importance of the 
modern X-ray monitoring missions like MAXI and the upcoming 
WFM onboard LOFT. 
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